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ABSTRACT 

We present polarimetric differential imaging (PDI) data of the circumstellar disk around the Herbig 
Ae/Be star HD100546 obtained with VLT/NACO. We resolve the disk in polarized light in the H 
and Ks filter between ~0.1-1.4" (i.e., ^10-140 AU). The innermost disk regions are directly imaged 
for the first time and the mean apparent disk inclination and position angle are derived. The surface 
brightness along the disk major axis drops off roughly with S{r) oc but has a maximum around 
0.15" suggesting a marginal detection of the main disk inner rim at ~^15 AU. We find a significant 
brightness asymmetry along the disk minor axis in both filters with the far side of the disk appearing 
brighter than the near side. This enhanced backward scattering and a low total polarization degree 
of the scattered disk flux of 14lg^% suggests that the dust grains on the disk surface are larger than 
typical ISM grains. Empirical scattering functions reveal the backward scattering peak at the largest 
scattering angles and a second maximum for the smallest scattering angles. This indicates a second 
dust grain population preferably forward scattering and smaller in size. It shows that, relatively, in 
the inner disk regions (40-50 AU) a higher fraction of larger grains is found compared to the outer disk 
regions (100-110 AU). Finally, our images reveal distinct substructures between 25-35 AU physical 
separation from the star and we discuss the possible origin for the two features in the context of 
ongoing planet formation. 

Subject headings: stars: pre-main sequence, stars: formation, planetary systems: protoplanetary disks, 
individual objects: HD100546 
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1. INTRODUCTION 

The direct detection and characterization of exoplan- 
ets around young A- type stars (e.g., Marois et al. 2008, 
2010; Kalas et al. 2008; Lagrange et al. 2009; Lagrange 
et al. 2010; Hinz et al. 2010; Quanz et al. 2010) have 
demonstrated that planets with several Jupiter masses 
can form in a broad range of orbital separations provid- 
ing important constraints on planet formation processes. 
Understanding the formation of planets in more detail re- 
quires a solid knowledge of the geometrical, physical and 
chemical parameters of the circumstellar disks in which 
planets are thought to form (e.g., Alibert et al. 2011). 

To study these planet-forming disks by means of direct 
imaging in the near-infrared requires high contrast tech- 
niques that sufficiently suppress the stellar light and, in 
case of ground-based observations, minimize the resid- 
ual speckle noise caused by uncorrected aberrations in 
the wavefront. A powerful technique to achieve this is 
polarimetric differential imaging (PDI; see, e.g., Kuhn 
et al. 2001; Potter 2003; Apai et al. 2004; Hinkley et al. 
2009) that takes advantage of the fact that the direct stel- 
lar light is essentially unpolarized whereas scattered light 
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from dust grains in the surface layer of a circumstellar 
disk is polarized (e.g., Fischer & Henning 1995; Fischer 
et al. 1996). In the past, this technique has been suc- 
cessfully applied to study the circumstellar disks around 
young massive stars (Jiang et al. 2005, 2008), interme- 
diate mass Herbig Ae/Be stars (e.g., Hashimoto et al. 
2011; Perrin et al. 2004; Oppenheimer et al. 2008; Hales 
et al. 2006), the TTauri star TW Hya (Apai et al. 2004; 
Hales et al. 2006) and some debris disks (e.g., Hinkley 
et al. 2009). 

Polarimctry from space with the Hubble Space Tele- 
scope (HST) has the advantage of a more stable point 
spread function (PSF) and a higher surface brightness 
sensitivity compared to ground-based facilities and the 
disk structure and dust grain properties have been in- 
vestigated for some selected objects (e.g., AU Mic and 
AB Aur; Graham et al. 2007; Perrin et al. 2009). Al- 
though most of the studies mentioned above had to cope 
with technical or observational limitations, such as the 
use of a coronagraph (which often reduces significantly 
the accessible inner working angles around the star), a 
rather modest telescope size or only single filter images, 
they provided valuable scientific insights for the observed 
targets. 

In this paper, we present the first PDI study of the 
transition disk around the well-studied Herbig Be star 
HD100546 (see, Table 1 for stellar parameters). The cir- 
cumstellar environment of this object has been inten- 
sively studied by previous imaging campaigns using di- 
rect imaging from the ground (Pantin et al. 2000, J and 
K band) and from space with HST NICMOS {l.Q^lm), 
STIS, and ACS (F435W, F606W, F814W) (Augereau 
et al. 2001; Grady et al. 2001; Ardila et al. 2007). In 
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particular, the scattered light images across optical and 
NIR wavelengths of the latter studies revealed a com- 
plex circumstellar structure interpreted as a large disk 
{R > 300 AU) and a remnant envelope. Emission at 
lO/xm possibly probing the warmer dust layers of the cir- 
cumstellar disk has been spatially resolved on tens of 
AU scales by means of mid-infrared (MIR) interferom- 
etry (Leinert et al. 2004; Liu et al. 2003). MIR spec- 
troscopy studying the emission features of heated dust 
grains on the disk surface hinted towards noticeable dust 
processing, i.e., a significant amount of crystalline grains 
and comparatively large grains (van Boekel et al. 2005). 
In particular, the resemblance of the MIR spectrum to 
that of the comet Hale Bopp (Malfait et al. 1998) led 
to the conclusion that dust processing similar to that in 
our own solar system must have taken place (Bouwman 
et al. 2003). Furthermore, to account for a deficit in the 
NIR excess compared to the strong MIR excess emis- 
sion, Bouwman et al. (2003) suggested the presence of 
an inner cavity with a radius of 10 AU in order to fit the 
spectral energy distribution. These authors speculated 
that a giant planet orbiting at roughly 10 AU may have 
opened such a gap. An inner cavity was further sup- 
ported by FUV long slit spectroscopy using HST/STIS 
(Grady et al. 2005) and spatially resolved echelle spec- 
troscopy of the [01] 630 nm emission line (Acke & van den 
Ancker 2006). The latter authors argued that temporal 
changes in their [01] line profile may be another signpost 
for the existence of the yet unseen massive planet orbit- 
ing within the cavity. From rovibrational CO emission 
lines Brittain et al. (2009) found further evidence for an 
inner cavity (^^13 AU radius) existing not only in the 
dust but also in the gaseous component of the disk (at 
least in CO). Qualitatively the same was found by van 
der Plas et al. (2009) although these authors could only 
put a lower limit of 8 AU on the inner radius of the CO 
emission region. Panic et al. (2010) used CO emission 
lines in the (sub-)mm to constrain the mass, kinemat- 
ics and outer radius of the molecular gas and found that 
more than 10~^ Mq of molecular gas is present in a Ke- 
plerian disk with an outer radius of ~400 AU. A wealth 
of additional molecular lines (e.g., CO, OH) as well as a 
69 fim silicate feature was found in data taken with the 
PACS instrument onboard the Herschel Space Observa- 
tory (Sturm et al. 2010). 

Currently, the most sophisticated model of the 
HD100546 system was presented by Benisty et al. (2010) 
who inferred the existence of an additional inner disk 
from _ftr-band interferometry data. These authors model 
the SED and the NIR/MIR visibilities with a tenuous in- 
ner disk component (reaching from ~0.3-4 AU) consist- 
ing of micron-sized grains, a gap (stretching from ^4 13 
AU) and a massive optically thick main disk (^13 350 
AU) . A sketch of the basic disk geometry is shown in Fig- 
ure 1. In this model, most of the disk mass is contained 
in the main disk comprised of mostly large dust grains (1 
/im to 10 mm), but a population of small grains (0.05 to 
1 /xm) is also required on the disk surface at a few tens 
of AU to account for the MIR and FIR excess emission. 

Using PDI on an 8m class, AO-assisted telescope in 
two infrared filters {H and Kg) without any occulting 
mask, our images allow us to analyze disk regions with 
unprecedented spatial resolution and small inner work- 
ing angle (~0.1" to the central star) which are difiicult 



to probe with other direct imaging techniques. Not only 
are we able to probe the disk geometry and structure but 
we can also put new constraints on the dust grain proper- 
ties on the disk surface between ^15-140 AU. In section 
2, we describe the observations. The data reduction is 
explained in detail in section 3. The key results are pre- 
sented in section 4 and analyzed in a broader context in 
section 5. Finally, we discuss our findings in section 6 
and summarize and conclude in section 7. 

2. OBSERVATIONS 

The observations were carried out on April 7, 2006, 
with the AO-assisted, high resolution NIR camera NACO 
(Lenzen et al. 2003; Rousset et al. 2003) mounted on 
UT4 at the VLT. The observing conditions were photo- 
metric. We used the SL27 camera with a pixel scale of 
27 mas/pixel and the detector was set to HighDynamic 
mode and read out in Double RdRstRd mode. NACO 
is equipped with a polarimetric focal plane mask, a ro- 
tatable half-wave retarder plate, and a Wollaston prism 
which can be inserted into the light path for imaging 
polarimetry. The Wollaston splits the I± and /|| polar- 
ization into the ordinary and extraordinary beam off- 
set by 3.5" in y-direction on the detector. Overlap 
of the two images is avoided by the polarimetric mask 
which restricts the original field of view of 27"x27" 
to X = 27" long and y = 3.1" wide stripes separated 
by 3.5" in r/-direction. For linear polarization measure- 
ments different polarization directions can be selected 
with the rotatable half-wave plate. One full polariza- 
tion cycle includes observations in four different retarder 
plate positions (0.0°, -22.5°, -45.0°, -67.5°). To correct 
for bad pixels we moved the object along the detector's 
X-axis between consecutive polarization cycles. In to- 
tal, we observed HD100546 in four different filters {H, 
Ks, NBl.QA, /B2.18). While we saturated the core of 
the PSF in the H and Ks filter so that the innermost 
5-6 pixels (in diameter) were no longer in the linear de- 
tector regime (i.e. > 10000 counts), the images in the 
narrow band filters were never saturated and only indi- 
vidual frames showed count rates slightly above the lin- 
earity threshold. The FWHM of the PSF in the narrow 
band filters was typically 4.0-4.5 pixels wide. 

In addition to HD100546, we also observed the main 
sequence star HD106797 as a reference star in the two 
broadband filters to check the validity of our data reduc- 
tion process. The same star was already used as a PSF 
reference star for HD100546 by Augereau et al. (2001). In 
the meantime, however, Fujiwara et al. (2009) detected a 
hot debris disk around HD 106 79 7 based on mid-infrared 
observations with the AKARI satellite making this ob- 
ject no longer an ideal reference source and we would 
pick a different object next time. 

While Table 1 summarizes the basic parameters of the 
two objects Table 2 provides a detailed summary of the 
observations and the observing conditions. 

3. DATA REDUCTION AND CALIBRATION 

3.1. Basic da,ta reduction steps 

The data reduction was carried out using self-written 
IDL- scripts unless otherwise described. Each exposure 
was first dark and flatfield corrected and then cleaned 
for bad pixels using a 3-ct filtering process that replaces 
the central pixel of a 5 x 5 pixel box with the mean of the 
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other pixels. Dither positions where the AO correction 
was poor (e.g., due to an open loop) were disregarded 
from any further analysis (see Table 2). Afterwards, the 
ordinary and extraordinary image in each exposure was 
cut out and all images of each filter (i.e., all dither po- 
sitions, all rotation angles) were aligned to one reference 
images that was chosen to be the first image of the sc- 
ries. As pointed out by previous authors (e.g., Apai et al. 
2004), the alignment of the images has to be done with 
sub-pixel accuracy in order to avoid the introduction of 
any spurious features and to eff'ectively suppress the un- 
polarizcd light from the central star. For the alignment, 
we used the IDL-based program IDP3^ that was initially 
developed for HST imaging data (Schneider & Stobie 
2002). To determine the centroid of each individual im- 
age, a two-dimensional Gaussian was fitted to the PSF 
over a 20-pixel radiiis. The shifting onto the reference 
image was carried out using "bicubic sine" interpolation. 
Once all images were aligned, we flagged those regions in 
the individual images where the count rates exceeded the 
linear regime, i.e., 10000 counts in this detector read-out 
mode. This way, in the final images at the end of the 
data reduction process non-linear pixels are flagged and 
disregarded in the analyses. 

3.2. Computation of Stokes parameter 

We measured the fractional Stokes parameters pq = 
Q/I and pu = U/I using double ratios as described in 
detail in Tinbergen (1996) or as applied in the polari- 
metric analysis of planets by Schmid ct al. (2006). 

From the two exposures taken with rctarder position 
0° and —45° we derive pq. This rotation of the rctarder 
plate allows to swap the two polarization images in the 
two WoUaston beams so that diff'erential effects cancel 
out. First, we derive 
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Here, the superscript denotes the rotation angle of the 
retarder plate, while the subscript denotes either the or- 
dinary or the extraordinary beam. By first dividing the 
ordinary images by the extraordinary images for a given 
retarder position, all time dependent factors such as at- 
mospheric transmission and seeing are cancelled out as 
these images were taken simultaneously. Then, by divid- 
ing the ratios of the different rctarder positions, all polar- 
ization influencing factors such as reflection coefficients 
of mirrors or individual pixel sensitivities are cancelled 
out. The fractional polarization Q/I is derived via 
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Q Rq 

^^2 = 7 = R^ 



(2) 



Since Q denotes the flux measured for the Stokes Q pa- 
rameter and / refers to the total intensity of the object 
(i.e., all counts summed up: /ord + 4xtra)- Pq can be 
expressed in \mits of "percent" or simply as a fraction. 
Using the images taken at the other retarder plate po- 
sitions (-22.5° and -67.5°) the fractional polarization for 
the U component, pu, is derived similarly. We computed 
Pq and pu for each dither position individually. After- 
wards the resulting images are averaged to increase the 

^ http://mips.as.arizona.edu/MIPS/IDP3/ 



signal-to- noise ratio (S/N) of the final image. The total 
fractional polarization pi is derived from the final pq and 
Pu images via 

' (3) 



PI 



To compute the individual Stokes components Q and 
U for each dither position we also applied the double 
difference approach as presented in Hinkley et al. (2009): 
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Q=(+Q-(-Q))/2. 
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Again, the superscript denotes the rotation angle of the 
retarder plate, while the subscript denotes either the or- 
dinary or the extraordinary beam. We refer the reader to 
the aforementioned paper for a detailed discussion of this 
method and its advantages. The Stokes U component is 
obtained in the same way using the images with rotation 
angles of -22.5° and -67.5°. Finally, the images from each 
dither position were averaged to increase the S /N of the 
final Q and U images. The intensity of the polarization 
(i.e., the total polarized fiux) is obtained from the final 
images via 

P = ^Q^ + m (7) 

3.3. Correction for instrumental polarization 

Being a Nasmyth instrument and never optimally de- 
signed for polarimctry, NACO suffers from significant 
instrumental polarization effects. Unfortunately, the in- 
strumental polarization changes whenever the positions 
of reflecting surfaces in the light path change relative to 
each other (i.e., due to different telescope pointing, differ- 
ent rotation angle of the camera, etc.). These effects have 
been investigated by Witzel et al. (2011) for the ii'-band. 
According to this study the introduced polarization is 
about 0.9 % for the telescope and 0.6 % for NAOS, while 
the C/ — > y cross talk between linear and circular polar- 
ization in NAOS can reduce the (7-polarization flux by 
up to 15 %. In this sense the measured U fluxes should 
be regarded as lower limits. At the time of our obser- 
vations the polarimetric characteristics of VLT/NACO 
were not known. For this reason we have chosen a cal- 
ibration approach which uses the central source as zero 
point polarization reference to estimate the amount of 
instrumental polarization for each data set. A detailed 
description is provided in Appendix A. 

3.4. Alignm,ent and photom,etric calihraMon 

It turned out that for data taken before fall 2009 the 
retarder plate position values saved in the image head- 
ers were off by ~13.2° compared to the actual angular 
positions on the sky due to a bug in the encoder posi- 
tions (Witzel ct al. 2011). This angular offset had to be 
subtracted from the camera positions saved in the image 
headers to obtain the correct orientation of the measured 
polarization direction on the sky (see. Table 2). Since we 
chose to align the disk major axis as observed with HST 
(Augereau et al. 2001) with the x-axis of the detector'' 

^ Due to an error in the initial rotation command the camera 
was rotated by 90° in the H band compared to the other filters so 
that the disk minor axis was aligned along the detector's x-axis. 
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the positive Q component is not aligned in north-south 
direction in our final images as it is the usual convention. 

The photometric calibration was based on the unsatu- 
rated images in the NBl.GA and 7S2.18 filters because 
the images in the broadband filters were saturated as de- 
scribed above. For all dark, flatfield and bad pixel cor- 
rected images we cut out the ~19.8" x 3.1" sub- images 
of the detector that contained either the ordinary (/q^j) 
or the extraordinary image (/cxti-a) of the star. The size 
of the sub-images along the x-axis was chosen so that the 
edges of the detector, where dark and white fringes were 
apparent, were excluded (^--^3.5" on both sides). The size 
along the y-axis was determined by the size of the Wol- 
laston mask. We computed and subtracted the local sky 
background in the individual sub-images. The remaining 
count rates in the sub-images were then regarded as the 
photometric signal from HD100564. 

To estimate the total intensity of the object we 
summed up the contribution from ordinary and extraor- 
dinary sub-images, i.e., 1^°^ ^ 1^°^ + I^^^^^, for each 
retarder plate and dither position and computed the 
mean value. For this we found approximately (1.446 ± 
0.012) xlO^ counts and (8.210 ± 0.079) x 10*^ counts in the 
A''B1.64 and 752.18 filter, respectively. The errors are 
the standard deviation of the mean (^0.8% and ^-^1.0%, 
in the respective filters). Using the transmission curves 
of the filters* we calculated that the A'' 731.64 filter has 
a transmission of only ~3.4% relative to the H band fil- 
ter. Between the 7B2.18 and the Kg filter this fraction 
amounts to ^13.2%. Taking these factors into account 
as well as the different integration times between the fil- 
ters (Table 2) we could estimate the expected count rates 
for the broadband filters. This intrinsically assumes that 
HD 100546 has the same magnitude in the A'^Bl.64 and 
7B2.18 filter as in the H and Ks filter, respectively (see, 
Table 1). Knowing now the count rates in the H and 
filters and the respective magnitudes we can derive the 
zero point from 

m2MASS = -2.5 ■ log(^) (8) 

with 7cR being the count rate and 7zp being the zero 
point. With the zero points we can then convert the 
count rates in each pixel in the final images to magni- 
tudes using again equation (8) above, and, finally, we 
can convert the magnitudes into siirface brightness 

S = m2MASs + 2.5 • log(A) [mag arcsec"^] (9) 

where A equals the square of the pixel scale (given in 
arcsec pixel" ^) of the camera. 

Taking into account all the difi^erent sources of uncer- 
tainties and errors as described in this and the previous 
sub-section, some of which arc difficult to quantify ex- 
actly, we estimate that the absolute fiux calibration of 
the final polarized flux images P is only good to 30-40%. 
The main uncertainty arises from the correction of the 
instrumental polarization. 

4. RESULTS 

4.1. Detection of polarized flux from HD100546 

* http;//www.eso.org/sci/facilities/paranal/instruments/naco/ 
inst / filters .html 



In Figures 2 and 3 we present the final Stokes Q and U 
images and the final pQ and pjj images in the 77 and Kg 
filter, respectively. A "butterfiy" pattern^ is apparent in 
all the images which we interpret as polarized emission 
coming from light that was scattered on dust grains on 
the surface layer of the circumstellar disk surrounding 
HD100546. Such a pattern is a manifestation of tangen- 
tial polarization vectors and similar images have been 
obtained for other sources, see e.g., Apai et al. (2004) 
and Hales et al. (2006). In Appendix B we provide an 
estimate for the typical signal-to-noisc ratio (S/N) of our 
images showing that the observed featiucs arc real. 

In our data, the corresponding H and Ks are qualita- 
tively very similar and, as expected, the pattern in the 
U images is rotated by roughly 45° compared to the Q 
images. In Figure 4 we show the final Stokes Q images 
taken in the NBIM and 752.18 filter. Again, the ob- 
served butterfiy pattern is similar to those obtained with 
the broadband filters, but due to shorter total integra- 
tion times and a more narrow bandwidth of the filters 
the S/N in these images is lower. This explains why the 
polarization pattern can not be traced to the same an- 
gular separation from the central star as in Figures 2 
and 3. Since our main analyses are based on the higher 
S/N broadband images we do not show and analyze the 
complete set of the images obtained in the A'^Sl.64 and 
7_B2.18 filters and mention only that the observed struc- 
tures agree very well. 

Figure 5 is identical with Figure 3 but shows the results 
for the reference star HD106797. No strong and con- 
sistent polarization pattern is apparent as most promi- 
nently seen in the pq and pij images. As mentioned in 
section 2, Fujiwara ct al. (2009) found that HD106797 is 
surrounded by a hot debris disk. They estimated the in- 
ner radius to be ~14 AU and put an upper limit of 41 AU 
for the outer radius. Thus, theoretically, we could have 
resolved this disk with our observations as HD106797 
and HD100546 are at the same distance. However, the 
HD106797 debris disk has so far only been detected at 
thermal infrared wavelengths. No scattered light images 
in the optical or NIR exist. Thus, given what is known 
about this source, it is impossible to estimate by what 
factor the polarized flux potentially coming from the de- 
bris disk is below our detection limits. However, as there 
are a lot of examples where bright debris disks were not 
detected in scattered light even with HST (Krist et al. 
2010) and taking into account that only a fraction of the 
scattered light will be polarized, it seems seems not sur- 
prising that the HD106797 debris disk is not detected in 
our data. Thus, we are convinced that our data reduc- 
tion process works correct and does not introduce any 
spurious features (see also. Appendix A and B where 
we discuss the polarimctric calibration and the S/N es- 
timates for both objects). 

A first more quantitative analysis of the HD100546 
butterfiy pattern is shown in Figure 6 where we plot the 
azimuthal profile of the polarization pattern in Stokes 
Q and pQ for the H band as well as Stokes U and pu 
for the Kg band. For four difii'erent annuli centered on 

^ The term 'butterfly" pattern arises from the fact that we see 
two bright lobes along one axis superimposed on two dark lobes 
rotated by ~90°. From equations (4)-(6) it is clear that both the 
bright and the dark lobes contain a polarization signal because we 
are looking at a difference image. 
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the central star the mean count rate was computed as a 
function of position angle (averaged over 10° bins). The 
error bars are the standard deviation of the mean values 
computed in the individual images at each dither posi- 
tion divided by the square root of the number of images 
that were combined (see, Table 2). 

The strongest sinusoidal variation is found between 
0.1-0.25" (corresponding to roughly 10-26 AU in pro- 
jected separation without taking into account disk incli- 
nation), but the pattern can be traced outward to angu- 
lar separations ~1.4" as can be directly seen in the pg 
and pu images of Figures 2 and 3. The polarization pat- 
tern in our two filters is qualitatively very similar, but 
quantitative differences between the filters become ap- 
parent. These differences will be discussed in more detail 
below where the flux calibrated images are analyzed. 

4.2. Disk orientation 

To assess the disk orientation we fitted isophotal el- 
lipses to the photometrically calibrated total polarized 
flux images (P) obtained in the H and Kg filter (upper 
panels Figure 7). The fit yields the apparent disk posi- 
tion angle and its inclination assuming the disk is intrin- 
sically centro-symmetric and geometrically thin (i.e., no 
disk flaring). For the fit we used three surface bright- 
ness ranges in the inner disk region (< 0.5") going from 
9.0-9.5, 9.75-10.25, and 10.5-11.0 mag/arcsec^, respec- 
tively. Fitting three different brightness ranges allows us 
to derive a more robust estimates about the disk orien- 
tation. For the disk position angle the fit yielded a mean 
position angle of 136.0°±2.9° (east of north) in the H 
band and 140.0°±2.4° in the Kg band. For the inclina- 
tion the values were 45.9°±1.7° in H and 48.0°±2.0° in 
Kg- The errors are the standard deviation of the values 
measiired in the three regions. Thus, the apparent disk 
orientation on these scales is comparable between the two 
filters and for the rest of the paper we adopt mean values 
of 138.0°±3.9° and 47.0°±2.7° for the position angle and 
the disk inclination, respectively. We note, however, that 
the true geometrical disk inclination might be somewhat 
different from this value as we used regions of the same 
brightness in polarized flux for the fit. These regions do 
not necessarily correspond to identical geometric regions 
on the disk surface as the polarization depends on the 
scattering properties of the underlying dust grains and 
the true disk flaring angle. However, most previous stud- 
ies used isophotal fitting to derive estimates about the 
disk orientation and we compare our values to previous 
measurements in Table 3. Our results agree with most of 
the previously derived values although all of them were 
obtained at larger angular separations. It is interesting 
to note that we obtain significantly different values if we 
fit isophots between ~0. 7-1.0" (13.0-13.5 mag/arcsec^). 
In this case the inclination is only ~ 20° and ~ 23° in 
H and Kg, respectively. Thus, at these separations the 
disk appears much more face-on but we note that the ob- 
served polarized fiux is no longer symmetric with respect 
to the disk major axis obtained from the inner disk (see. 
Figure 7) and a more detailed modeling of the scattering 
properties of the dust grains is required. 

4.3. Disk surface brightness 
4.3.1. Radial profiles and asymmetries 



We extracted the profiles of the polarized fiux P and 
the polarization fraction pi along the semi-major and 
semi-minor axes of the disk in both filters (see. Figure 7). 
For both filters we used a 3-pixel wide slit with an ori- 
entation of 138° for the major axis (Table 3) and a per- 
pendicular orientation for the minor axes. The resulting 
mean profiles are shown in Figure 8. The error bars are 
the standard deviation in the 3-pixel wide slit for a given 
position. The disk has an overall higher polarized fiux 
surface brightness S{r) in the Kg band than in the H 
band. However, in both filters the brightness decreases 
similarly fast as a function of radius. Between 0.1-1.4" 
we find the following power-law exponents that fit the 
profiles best along the major axis: —3.47 ±0.05 (H) and 
—2.87 ± 0.04 (Kg) along the SE semi-major axis, and 
-2.84 ± 0.04 (H) and -2.86 ± 0.01 (A',) along the NW 
semi-major axis. The fits were " error- weighted" to ac- 
count for the uncertainties in the measured radial profile. 
The errors quoted for the power-law indices are fitting 
errors. With one exception the profiles show the same 
radial dependency, roughly S{r) oc r~^'^^. Only in the 
south-east the flux drops even steeper in the il-band with 
Sir) oc r-3-47_ 

In case of a centro-symmetric disk seen face-on the 
brightness proflle looks identical in all directions. For 
an inclined disk we would expect the radial proflles to 
be identical on both sides of the star along the disk ma- 
jor axis (as it is observed in the Kg filter). Along the 
disk minor axis, which is an axis of mirror symmetry 
in case of an inclined but intrinsically centro-symmetric 
disk, the forward and backward scattering properties of 
the dust grains are directly imprinted and the profiles 
will be different for both semi-minor axes. This effect 
is seen in Figure 7: The polarized fiux drops much more 
rapidly between roughly 0.1-0.8" in the south-west direc- 
tion compared to the north-east direction. This trend is 
observed in both filters. It is even more easily seen in the 
brightness profiles for the polarization fraction where the 
fractional polarization in the south-west is significantly 
lower than that in the north-east (bottom panels of Fig- 
ure 7). 

4.3.2. Disk polarization fraction 

The fractional polarization pi measured in our data 
does not measure the fractional polarization of the scat- 
tered light coming from the disk alone as it relates the 
observed polarized fiux to the total intensity composed 
of both the flux coming from the disk and that coming 
from the PSF of the central star. However, by combining 
the results from Augereau et al. (2001) with ours we can 
get an idea of the polarization fraction of the light scat- 
tered from the disk surface. Augereau et al. (2001) found 
an azimuthally averaged value of roughly 10.9±0.5 mag 
arcsec"^ for the surface brightness of the disk in total 
scattered hght using HST/NICMOS at 1.6 fim. Aver- 
aging the polarized flux along the semi-major and semi- 
minor axes at the same separation in our _ff-filter image 
with comparable central wavelength (Figure 8), we find a 
value of roughly 13.0±0.4 mag arcsec^^. This leads to a 
total polarization fraction of the disk of ~14t8^% where 
the errors relate to the minimum and maximum values 
given the photometric uncertainties. Ideally, to reduce 
systematic errors, such a computation is based on data 
where the total disk intensity /oisk and the polarized flux 



6 



Quanz et al. 



P is measured at the same time as it is possible with HST 
(e.g. Perrin et al. 2009). 

4.4. Disk color 

Having the surface brightness of the polarized flux 
available in two filters we can compute the [H — K s\-co\ov 
of the disk surface as shown in Figure 9. The image has 
been smoothed with a moving box of 4x4 pixels so that 
large scale structures are enhanced. Filtering the image 
with a 2-D Gaussian kernel with a FWHM of 4 pixels in- 
stead of the moving box does not change the results and 
leads to an almost perfectly identical image. It shows 
that the disk surface has on average a red color where 
the reddest parts with [H — Kg] > 0.5 mag arcsec"^ 
run roughly along the disk major axis. The innermost 
1.0" (in diameter) of this area appear like a bar with 
still moderate colors {[H — Ks] = 0.5...1 mag arcsec"^), 
while fmther out the disks becomes increasingly redder 
over a bigger wedge of the disk surface. It is interest- 
ing to note that the position angle of this red symme- 
try axis is ^^160° (east of north) which is ^-^20° larger 
than the position angle derived from the isophot fitting 
in section 4.3.1 but consistent with the disk position an- 
gle found by Augereau et al. (2001) in total intensity at 
1.6 fxin. Perpendicular to this axis, i.e., roughly along 
the disk minor axis of the disk, the color is typically 
[H — Kg] < 0.5 mag arcsec"^. In this context, we remind 
the reader that the intrinsic [H — Kg] color of a B9V star 
is ~0 mag and that the NIR excess emission and the ob- 
served NIR colors of IID100546 (see Table 1) have their 
origin most likely in the tenuous inner disk of the system; 
Benisty et al. (2010) found from their model that only 
'~20% of the total observed if -band flux comes directly 
from the star, while direct thermal emission from the in- 
ner disk and thermal emission that was scattered within 
the inner disk account for ^38% each. This leaves ~4% 
of the total flux to be coming from the star and then 
being scattered from the inner disk. Hence, most of the 
scattered (polarized) flux that we observe here was emit- 
ted and/or scattered from the inner disk. 

4.5. Sub-structure in the innermost regions 
4.5.1. A hole and a clump? 

In Figures 10, 11 and 12 we analyze the structure of 
the disk in more detail. These images have been rotated 
clockwise so that the disk major axis runs horizontally. 
Image 10 shows the high-pass filtered version of the P 
images shown in Figure 7. Here, we have subtracted a 
smoothed version (smoothing box size 10x10 pixels) of 
the images from the images themselves so that large scale 
structures are filtered out and small scale structures are 
enhanced. It shows that outward of ~0.5" (50 AU) no 
sub-structures are apparent and the disk surface appears 
rather homogeneous and smooth. In the inner parts of 
the disk, however, two features appear as indicated by 
the arrows and numbers in Figure 10: 

(1) Up-right of the image center, there seems to be a 
"hole" in the disk. To explore this region further, we 
zoom in the innermost ~40 AU of the disk and show the 
polarized flux images as well as the fractional polariza- 
tion images in the H and Kg filter in Figure 11. It shows 
that in the north-east of the disk (at a position angle of 
f«12° E of N) we have a signiflcant drop in polarized flux 



and fractional polarization which shows up most clearly 
in the fractional polarization images. The minimum of 
the polarized flux is located at the same location in all im- 
ages, i.e. at w8 pixels radial separation from the central 
star (corresponding to «22 AU projected separation). 
Also, the magnitude of the flux drop is similar in all im- 
ages. In case of a centro-symmetric disk, the left-hand 
and right-hand side of the images should be identical as 
the disk minor axis of the disk runs vertically through the 
center. Hence, we can use the mirror-symmetric location 
in the left-hand side of the images to quantify the flux 
drop on the right-hand side. It shows that we measure 
only 69% ± 4% of the flux / fractional polarization in an 
aperture with 3-pixels radius. As this figure is the mean 
and the standard deviation of the ratios in all 4 images 
shown in Figure 11 we conclude that this flux drop is real 
as it occurs at the same position in H and Kg- 

(2) Extending almost exactly to the north-west along 
the semi-major axis of the disk we find a local flux max- 
imum in the disk (a 'clump'). Although this feature is 
present in both filters the location of its peak flux is not 
exactly identical. The separation from the central star 
amounts to ^0.35" in both images but the position an- 
gle (E of N) is roughly -21° in the H filter and -32° 
in the Kg filter. Still, the region of enhanced polarized 
fiux clearly overlaps in H and Kg. In Figure 12 we have 
subtracted the mirrored left-hand side of the disk from 
the right-hand side, to check whether this apparent fiux 
maximum in Figure 10 is real. As mentioned above, an 
inclined centro-symmetric disk should have a mirror sym- 
metry with respect to its minor axis. However, we can 
clearly see that the right-hand side of the disk has a flux 
excess in the inner ^-^0.6" compared to the left-hand side 
in both filters, and the 'clump' seen in Figure 10 is part 
of a larger structure. At larger radial distances no ad- 
ditional structures are seen and both sides of the disk 
appear very similar and cancel out in the subtracted im- 
ages. 

Since the differential PDI technique compensates for 
instrumental effects within the CONICA camera and 
since we took our data with the camera rotated by 90° 
between the H and Kg filter we conclude that these fea- 
tures are real and cannot be a residual or ghost intro- 
duced by the camera, the telescope or the AO-system. 
Furthermore, we applied a second high-pass filtering 
method by taking the Fourier transform of the images 
and filtering all the low frequency components. Also here 
we recovered both features in both images. Finally, the 
features are located in image regions that show compar- 
atively high S/N ratios (S/N>10 in the individual Q and 
U images at the location of the features; see, Appendix 
B). We will discuss these features in more detail in sec- 
tion 6.3. 

4.5.2. The rim of the main disk 

As mentioned in the introduction, the current model of 
the HD100546 system locates the inner rim of the main 
disk at ~13 AU (Benisty et al. 2010) and there is ob- 
servational support for this from different studies (e.g., 
Grady et al. 2005). We checked in our images whether 
we are able to confirm this finding. Figure 13 shows a 
horizontal cut through the center of the images shown in 
the upper row of Figure 1 1 . This is essentially identical 
to the -0.4 to -1-0.4 regions of the plots shown in the up- 
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per row of Figure 8. We used a 3-pixcl wide horizontal 
box and computed the mean and the standard dciviation 
of the pixels for a given position. The innermost regions 
close to the star have been masked out if any of the pixels 
in the slit at this position showed count rates exceeding 
the linear detector regime. It shows that we find the 
maximum polarized flux between 0.1 0.15" on both sides 
of the central star in both filters and a steep drop of flux 
closer to the star. The angular separation of the max- 
imum flux corresponds to 10-15 AU at the distance of 
HD100546 which would support the idea that there is 
a jump in the density of scattering dust grains at this 
location. 

The interpretation of these results should be taken with 
caution. Although we find the maximum flux at the sus- 
pected location of the inner rim, there is only 1 radial 
pixel closer to the star on both sides and in both fil- 
ters that did not suffer from saturation effects. Also, 
only along the north-west axis (positive x-values) the 
flux drop appears to be signiflcant taking into account 
the error bars. However, the consistent behavior we find 
for the polarized flux profiles close to the star in both 
filters allows us to conclude that we may indeed have di- 
rectly imaged the suspected inner rim of the main disk 
between 10-20 AU. ideally, additional new data with a 
non- (or less) saturated PSF core would be required to 
fully confirm this interpretation. 

5. ANALYSIS 

5.1. Disk geometry - Backward scattering preferred 

In order to relate the observed polarized flux to some 
dust grain properties of the disk surface, it is necessary to 
understand which side of the inclined disk is the near side 
and which is the far side. Based on the rotation signa- 
ture of the [01] gas in the inner 100 AU, Acke & van den 
Ancker (2006) concluded that the disk is rotating coun- 
terclockwise, i.e., the north-west side is receding. This 
is also confirmed by the CO measurements from Panic 
et al. (2010). Combining this information with the large- 
scale spiral structure seen by Grady et al. (2001) and also 
Ardila et al. (2007), and assuming that these spiral struc- 
tures are trailing spiral arms, the north-east side of the 
disk is the far side and the south-west side of the disk 
the near side relative to Earth. Based on this, the ob- 
served brightness asymmetry along the disk minor axis 
(section 4.3.1) tells us that the far side of the disk ap- 
pears brighter in polarized flux than the near side. This 
in turn means that the dust grains responsible for the 
polarized flux are preferentially backward scattering. 

We note that the observed brightness asymmetry 
within the innermost «0.8" along the disk minor axis was 
already recognized in HST/NICMOS data by Augereau 
et al. (2001) and anisotropic scattering properties of the 
dust grains was mentioned as one possible explanation. 
Interestingly, in the HST/ACS data presented by Ardila 
et al. (2007) there is also a brightness asymmetry along 
the minor axis but in the opposite sense: the amount 
of scattered light in the south-west is larger than in the 
north-east in all thrcie optical bands. However, these 
authors were only sensitive to disk regions >1.6" as sub- 
traction residuals dominated in the innermost regions. 
Still, this result needs to be considered if one tries to 
derive dust grain properties from a combination of all 
scattered light images. We discuss this in more detail 



below in section 6.2. 

5.2. Polarized scattering functions 

To further quantify this brightness asymmetry we com- 
puted the polarized scattering function for the disk in 
both filters (Figure 14). For this purpose we used again 
the polarized flux images shown in the upper panels of 
Figure 7. Given the structures we found in the north- 
west side of the disk (sec, section 4.5), the scattering 
functions were computed for the opposite disk side. The 
scattering angle 9 \s a, function of the azimuthal angle 4> 
(starting from the semi-minor axis in the north-cast go- 
ing counterclockwise to the semi-minor axis in the south- 
west), the disk inclination i and the disk flaring angle 7 
and is given by 

6» = 90° + (i - 7) • cos (for(/-<90°) , . 
6* 90° + (i + 7) • cos (f) (for (j) > 90°) 

For the disk inclination wc adopted the mean value of 
47° (see Table 3) and for the flaring angle 7° which cor- 
responds roughly to the flaring index j3 used in the disk 
model from Benisty et al. (2010). For two different annuli 
(40-50 AU and 100-110 AU) we computed the mean and 
the standard deviation of polarized flux in wedges of 10° 
as a function of 9 and plotted the results in Figure 14. 
For the annuli we also had to take into account the disk 
inclination to make sure we were analyzing the same disk 
regions in terms of physical and not projected separation 
from the star. The curves were normalized to unity at 
the brightest disk region which is the far side of the disk 
as described above. The maximum scattering angle is 
roughly 130° while on the near side of the disk it is a bit 
less than 40°. 

In these scattering functions a few points points are 
noteworthy: (1) All four curves are qualitatively very 
similar: Going from the backward scattering peak at the 
largest scattering angles to smaller 9, all curves drop, 
show a minimum between 9 = 40°-70° (i.e., a scattering 
angle with minimum scattering efficiency) before they 
rise again for even smaller angles. (2) Comparing the 
curves for the H filter with those for the Kg filter, the 
flux drop seems less pronounced in the Kg fllter in both 
annuli. (3) Comparing the curves for the two annuli, in 
both fllters the black curves (probing 40-50 AU) seem to 
fall off less steep than the red curves (probing 100 110 
AU) going from the backward scattering peak to smaller 
values for 9. Also, the black curves fall more or less 
steadily imtil they reach their minimtun and the rise at 
the smallest scattering angles is not as high as for the 
red curves. In contrast, the red curves, after their sharp 
drop-off coming from the backward scattering peak, re- 
main almost flat until they increase again at angles <50°. 

6. DISCUSSION 

6.1. Brightness profile and comparison to other disks 

The measured profiles show that the scattered surface 
brightness along the semi-major axes in the sky plane 
drops off roughly S (x r^'^. If one includes disk fiaring 
then the expected profile becomes less steep and fits even 
better to our observed values. Taking into account the 
f(r) oc drop in stellar radiation this implies that 
the surface density of the scattering dust grains drops 
off roughly E oc r~^ if we assume constant scattering 
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properties of the dust grains over the observed radial 
range. 

Our values for the inner ~140 AU of the disk are 
in good agreement with the results from Augereau 
et al. (2001) who found an azimuthally averaged sur- 
face brightness profile S{r) oc r~^-® between ^55-250 
AU in their HST/NICMOS l.e^im data. Further out, 
the power law seems to become steeper: Augereau et al. 

(2001) measured S{r) oc r-^-^ between 250-350 AU while 
Ardila et al. (2007) found S{r) oc r'^-^ between 160-500 
AU in their optical HST/ACS data. This suggests ei- 
ther a discontinuity in the surface density of the dust 
disk, different dust grain properties at larger disk radii 
or a different disk geometry (e.g., changing flaring angle). 
We come back to this point below in the next section. 

Comparing the observed surface brightness proflle to 
other directly imaged objects we find that also the Her- 
big Ae star HD169142 shows S(r) oc r^^ between 80- 
190 AU (Grady et al. 2007, ; HST/NICMOS). For lower 
mass TTauri stars Apai et al. (2004) found S{r) oc r~^-^ 
for TW Hya between 50-80 AU using NACO/PDlio 
and Schneider et al. (2003) found S{r) oc j^-S ' -S-S f^j. 
GM Aurigae for distances >85 AU using HST/NICMOS 
data. Irrespective of stellar mass, many circumstellar 
disks show a very similar surface brightness profile in 
scattered light images for disk radii >50 AU consistent 
with optically thick (and geometrically thin or only mod- 
erately flared) disks. 

Interesting is the comparison of HD100546 with TW 
Hya as both objects are transition disks and both were 
observed with NACO/PDI. For HD100546 the strongest 
polarimetric signal is found in the innermost regions of 
the disk (i.e., in the annulus from 0.1-0.25" correspond- 
ing to roiighly 10-26 AU in projected separation without 
taking into account disk inclination) while Apai et al. 
(2004) found the stronger signal not close in but further 
out in the TW Hya disk (0.76-1.03" corresponding to 
roughly 43-58 AU). Thus, apparently those regions in the 
TW Hya transition disk that are optically thin and de- 
pleted of dust reach further out compared to HD100546. 
And since TW Hya lies at only ~54 pc distance, the 
data from Apai et al. (2004) seem to resolve these re- 
gions of reduced dust surface density. As we have seen 
in section 4.5.2 our data tentatively resolve the inner rim 
of the HD100546 main disk but probing even further in 
(i.e, closer than 10 AU) is beyond the capabilities of our 
observations. 

6.2. Dust grain properties 

Having disk images in polarized flux one can constrain 
dust grain properties much better than based on disk 
intensity image alone (e.g., Graham et al. 2007; Perrin 
et al. 2009). However, in order to interpret the polarized 
scattering functions shown in Figure 14 one has to keep 
in mind that they are the product of two functions which 
depend both on diist grain properties such as grain sizes, 
composition and porosity: The " scattering function" de- 
scribing the dust scattering efficiency as a function of 
scattering angle, and the " polarization function" describ- 
ing the introduced scattering polarization as a function 

1" Wc note that based on HST/NICMOS data Weinberger et al. 

(2002) found a more shallow profile for TW Hya of S{r) oc r-2-6 
between 45-150 AU. 



of scattering angle. 

As mentioned above and shown in Figure 14, we find 
that the far side of the disk appears brighter than the 
near side indicating that backward scattering if preferred 
over forward scattering. This, however, requires scatter- 
ing particles that are typically at least as large as the 
wavelength of the incoming light. The strong flux de- 
tected in mm-observations (Henning et al. 1998; Wilner 
et al. 2003) already indicated the existence of dust grains 
that are considerably larger than those found typically 
in the ISM (i.e., sub-micron sized). While most flux ob- 
served in the mm-domain is in general attributed to the 
mid-plane of a circumstellar disk, our scattered light ob- 
servations probe grain properties of the disk surface. 

To our knowledge there are only two other prominent 
examples of disks where backward scattering dominates: 
The first one is Sattirn's Rings which arc made up of 
icy particles that arc predominantly cm-sized or larger 
(see, Dones et al. 1993, and references therein). The sec- 
ond one is the debris ring around Fomalhaut. Min et al. 
(2010) showed that the scattering function observed by 
Kalas et al. (2005) could be matched with grains several 
tens of micrometer in size. Comparing our scattering 
functions to those published by Min et al. (2010) it seems 
as if grains of similar size could match our observations. 
With this interpretation one has to be cautious, though: 
First, the HD100546 is very likely optically thick which 
means that optical depths effects and multiple scattering 
will influence the scattering function compared to the 
optically thin, single scattering case. Secondly, we only 
observe the polarized component of the scattered light. It 
is important to consider that the polarization efficiency 
typically peaks for scattering angles close to 90° and is 
always zero in exact forward or backward scattering di- 
rection. Thus, a direct quantitative comparison with the 
dust grains sizes (Min et al. 2010) and dust composition 
of other circumstellar disks (Graham et al. 2007; Perrin 
et al. 2009) requires detailed modeling which is beyond 
the scope of the present paper. Qualitatively, however, it 
is clear that rather large, micron-sized grains are required 
to explain the strong backscattering peak we observe. 

In addition, the rather low polarization degree of the 
scattered light of ^14% (measured at a separation of 
0.5"; see, 4.3.2) also supports the existence of large grains 
on the disk surface as the polarization efficiency decreases 
with increasing particle size (Graham et al. 2007). The 
value wc find for HD100546 is lower than that of AB 
Aur where Perrin et al. (2009) found roughly 40% of the 
scattered light to be polarized (on average) in an annulus 
of 0.7-1.5" (~100-215 AU). For the TW Hya disk. Hales 
et al. (2006) found the absolute polarization to increase 
from -10% at 0.5" to -30% at 1.2" in both J and H 
band. Further out, at roughly 2", the value increased 
up to —40% in H but faUs back to -10% in J what 
the authors attributed to a steeper surface brightness 
profile of polarized flux in J. For HD100546, we can 
expected that, since our surface brightness profile along 
the disk major axis and that of Augereau et al. (2001) are 
similar, the ratio of polarized light to total scattered light 
does not change significantly and hence the polarization 
degree between 0.5-1.4" should not change significantly 
either. Perrin et al. (2009) used a disk model to fit their 
observations and the best fit model uses the typical grain 
size distribution with a power low index of -3.5 and a 
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maximum grain size of 1/im. The same maximum grain 
size is currently used in the latest disk models for the 
surface layer of the HD100546 main disk (Benisty et al. 
2010) but our data suggests that a significant fraction of 
large grains is required. 

Further, we can speculate that the differences we see 
in the scattering function between the inner disk and 
the outer disk (sec, Figure 14) could hint toward slightly 
different dust grain populations. Larger grains tend to 
have their scattering minimum at smaller scattering an- 
gles (Min et al. 2010). Perhaps the grains in the inner 
disk regions are slightly larger than those in the outer 
disk regions. This hypothesis is hirthcr siipported by 
the existence of the second scattering peak we observe 
for the smallest scattering angles: This peak is indicative 
of an additional dust grain population of smaller grains 
that are preferably forward scattering. Figure 14 shows 
that this peak is stronger for the outer disk regions (100- 
110 AU) compared to the inner disk regions (40-50 AU). 
In this context it is interesting to note that Augereau 
et al. (2001) derived a blow-out size for the dust grains 
of '^IS /im, which would mean that smaller dust grains 
should be blown out to larger radii and eventually out of 
the system. As these computations took only radiation 
and gravitational forces into account it is unclear how 
the total gas-to-dust ratio on the disk surface is and how 
potential pressure for the remaining gas (see, section 1) 
affects this number. However, having a higher fraction of 
small grains further out on the scattering surface could 
also explain the brightness asymmetry that Ardila et al. 
(2007) found along the disk minor axis in their optical 
HST/ACS images. In their data forward scattering, pos- 
sibly coming from smaller particles, was enhanced. This 
brightness asymmetry was most prominent for disk radii 
>2.0" (>200 AU). Inside of 1.6" the authors were limited 
by the HST/ACS coronagraph and subtraction residuals 
and no information could be derived. Thus, our data are 
not probing the same disk region and maybe different 
dust grain populations are involved: smaller dust grains 
(blown out?) at large separations, and an increasing frac- 
tion of larger particles closer to the star. We note that, 
similarly, another abrupt change in the color and bright- 
ness in their scattered light images for radial distances 
> 8" led Ardila et al. (2007) to the conclusion that they 
were probing yet again a different regime of the circum- 
stellar matter - the authors referred to it as an extended 
envelope - possibly featuring again a different dust grain 
population. 

6.3. Disk sub- structures as signposts of ongoing planet 
formation? 

The small inner working angle of our observations re- 
vealed structures in the inner disk regions of IID100546 
that have not been detected previously (see 4.5). Similar 
structures have recently been observed in the innermost 
regions of the AB Aurigae disk also using ground-based 
polarimetric imaging (Hashimoto et al. 2011). It is diffi- 
cult to identify the underlying physical reasons for these 
structures as different mechanisms can Iciave an observa- 
tional imprint on the disk surface density ranging from 
disk intrinsic phenomena such as gravitational instabili- 
ties or magneto-rotational instabilities to interactions be- 
tween the disk and an embedded body such as a young 
planet. 



Both features, i.e., 'the hole' and 'the clump', were 
detected in both filters in the p/ and P images so 
that we can be sure that they are features of the disk 
(see 4.5). However, while interpreting these features we 
have to keep in mind that our observations are prob- 
ing the dust component of the disk surface which is 
presumably optically thick at NIR wavelengths. We 
can not say whether the features are caused by a defi- 
ciency/overabundance of dust grains at these locations 
(i.e., a physical 'hole'/'clump' in the disk) or whether 
different scattering properties of the dust grains are re- 
sponsible for the observed features. Since both features 
are very localized and no additional large scale disk struc- 
tures are present in our data^^ we will discuss 'the hole' 
and 'the clump' in the context of interactions with em- 
bedded bodies and will not consider possible mechanisms 
driven by the disk alone or very localized changes in the 
grains' scattering properties. 

Numerical simulations have demonstrated that embed- 
ded planets can rapidly create gaps and sub-structures in 
the dust component of a circumstellar disk which depend 
on the mass of the planet but also on the size of the dust 
grains with more massive planets creating deeper and 
wider gaps and larger grains being more affected (e.g., 
Fouchet et al. 2010; Paardekooper & Mellema 2006). 
These features are most prominent in the disk midplanc, 
i.e., where the planet is formed, and are most easily seen 
in (sub-mm) observations. This in turn means, that a 
lower mass body orbiting in a disk is causing only lit- 
tle signatures and is more challenging to see in scattered 
light images (see, e.g.. Wolf 2008). However, the grav- 
itational influence of the planet could remove dust and 
gas particles locally and create small scale depressions 
on the disk surface. This could create planet-induced 
disk 'shadows' in the NIR which have been modeled by 
Jang-Condell (2009). Even if the parameter space stud- 
ied in the aforementioned paper does not cover the setup 
we observe here, it provides at least qualitatively an in- 
teresting explanation for 'the hole': Even bodies with 
masses of lOM® < M < 5OM0 can leave an observa- 
tional imprint on the disk surface on the order of the 
flux drop we observe here (Jang-Condell 2009). 

Since the apparent distance between 'the hole' and 'the 
clump' is rather small (^9 pixels) it seems interesting 
to investigate whether they may be somehow related. 
Taking into account the rotation sense of the disk, 'the 
clump' is trailing behind 'the hole'. If we de-project the 
apparent separations from central source (rdump ~ 13 px; 
'"hole ~ 8 px) considering the disk inclination angle and 
the position angle of the features we find that 'the clump' 
has a physical separation of ~35 AU while for 'the hole' 
the separation amounts to ^27 AU. In the above men- 
tioned numerical simulations from Fouchet et al. (2010) 
it is shown that on the outer edge of the planets or- 
bital gaps the dust surface density is enhanced. Also, in 
the Jang-Condell (2009) simulations, in addition to the 
planet- induced disk shadow roughly at the planet's lo- 
cation, the disk will appear brighter radially behind the 
planet compared to the surroundings. It would be inter- 
esting to see whether the static radiative transfer models 
by Jang-Condell (2009), when combined with dynamical 

As mentioned above, there is evidence for larger scale struc- 
tures such as spiral arms at greater radial distances in HST images. 
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SPH simulations, would predict a brighter spot trailing 
behind the 'shadow'. This would niateh with our obser- 
vations and support a physical link between 'the hole' 
and 'the clump'. 

We emphasize that our observations and the struc- 
tures we see in our data provide no immediate proof 
that planet formation is taking place in the disk around 
HD100546. The observed features are, however, consis- 
tent with the predictions of certain theoretical models 
dealing with the observational signatures of planet-disk 
interactions. A good observational test for our structures 
will be to rc-obscrvc HD100546 in the future. If 'the 
hole' and 'the clump' are caused by (an) unseen plane- 
tary companion(s) then a Keplerian orbital motion with 
a baseline of 10 years will lead to a rotation of roughly 30° 
(assuming circular orbits), enough to be detected even if 
the estimated disk inclination is taken into account. 

6.4. The big picture - HD 100546 a remarkable 
transitional disk 

In combination with previous findings, the results pre- 
sented here emphasize that the HD100546 circumstcUar 
disk seems to be an ideal laboratory to study the forma- 
tion of planets: (1) The inner regions between 4-13 AU 
are mostly devoid of dust and gas (e.g., Bouwman et al. 
2003; Grady et al. 2005; Acke & van den Ancker 2006; 
van dcr Plas et al. 2009; Bcnisty et al. 2010) possibly 
cleared out by an unseen (planetary mass?) companion 
(e.g. Acke & van den Ancker 2006); (2) the MIR spec- 
trum is remarkable similar to that of comet Hale-Bopp 
(Malfait et al. 1998; Bouwman et al. 2003) indicating 
that some fraction of the HD100546 dust component is 
similar in composition to diist foimd in our Solar Sys- 
tem; and (3) our observations revealed structures in the 
outer parts of the disk (at ~30 AU separation) that could 
be the signpost of an/additional body/bodies forming at 
these locations. 

At the same time, HD100546 seems to be still accreting 
(Vieira et al. 1999) from its inner disk and the main disk 
contains a reservoir of molecular gas that is still massive 
enough (>10^^ Mq) to form additional planetary bodies 
(Panic et al. 2010). It seems somewhat surprising that 
HD100546 still has such a significant gas reservoir, as 
gaseous disks around early type stars typically disperse 
on timescales <3 Myr (Kennedy & Kenyon 2009). Even 
if recent estimates of the surface gravity of HD 100546 in- 
dicate an age between 5 -10 Myr (Guimaracs et al. 2006) 
instead of the more frequently \ised age of >10 Myr (van 
den Ancker et al. 1997), the measured disk mass is ex- 
ceptional. Finally, it should be kept in mind that at 
larger separations (>2") the disk becomes highly struc- 
tured again as seen in the scattered light images of HST 
(Grady et al. 2001; Augereau et al. 2001; Ardila et al. 
2007). Quillen (2006) suggested that the spiral-arm like 
features and brightness variations detected in these im- 
ages result from a warped and twisted disk geometry out- 
side of ^200 AU. An interesting mechanism to produce 
such a tilt in the disk would be precession induced by of a 
swarm of planetesimal-like bodies embedded in the disk 
at larger separations (50-200 AU) (Quillen 2006). The 
presence of such bodies would also support the existence 
of larger dust grains at these separations (e.g., created 
through collisions) as suggested by our data. 

Overall HD100546 may be a great showcase example 



of a young star where we can directly witness the forma- 
tion of a new planetary system with at least one planet 
already orbiting within the disk gap and additional plan- 
ets and planetesimals still forming in the disk in which 
they are embedded. 

7. SUMMARY AND CONCLUSIONS 

In this paper we presented polarimetric differential 
imaging (PDI) observations of the circumstellar disk 

around the Herbig Be star HD100546 using the AO as- 
sisted high-resolution camera NACO at the VLT. We 
detected and clearly resolved the disk from ^0.10-1.4" 
(corresponding to the inner 10-140 AU) in polarized flux 
in the H and Kg filter. The innermost disk regions <50 
AU are directly imaged for the first time emphasizing the 
power of PDI. Our main findings can be summarized as 
follows: 

• The disk inclination measured from isophot fit- 
ting in H and between 0.2-0.5" is 47.0°±2.7° 
in agreement with previous estimates from direct 
imaging. For the disk position angle we find 
138.0°±3.9° again in general agreement with pre- 
vious work. 

• Fits to the surface brightness S{r) along the semi- 
major axes between 0.1-1.4" yield power laws with 
S{r) oc r-~^-* - ~^-^ in both filters. Only in the 

H band the south-east semi-major axis shows a 
steeper profile with S{r) oc r~^'^. 

• The brightness profile of the disk in polarized flux 
is asymmetric along the disk minor axis. In both 
filters, the far side of the disk appears brighter than 
the near side indicating that the dust grains on the 
disk surface are preferably backward scattering. 

• The above mentioned brightness asymmetry in 
combination with a low polarization fraction of 
the scattered light of only 14tg®% indicates that 
the dust grains probed with our observations have 
probably sizes on the order of micrometer, i.e., sig- 
nificantly larger than the size of pristine ISM dust 
grains. Grains of this size on the disk surface chal- 
lenge current theoretical models of the HD100546 
disk that typically assume a population of smaller 
grains. 

• The polarized scattering functions we measure in H 
and Kg cover scattering angles 9 « 40° 130° (tak- 
ing into account the disk inclination and a theoreti- 
cally motivated flaring angle) and are qualitatively 
very similar in both filters. While the maximum 
is found for the the largest scattering angles (the 
manifestation of preferred backward scattering) the 
minimum typically lies between 50°-70°. For even 
smaller scattering angles the functions do rise again 
which could be caused by an additional population 
of smaller dust grains which scattered preferably in 
the forward direction. Comparing the exact shape 
of the scattering function (i.e., position of the min- 
imum, intensity of the forward scattering peak) for 
disk regions between 40-50 AU to that between 
100 -110 AU suggests that the relative fraction of 
small grains is higher in the outer disk regions. 



NACO PDI observations of the disk around HD100546 
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• The small inner working angles of our imaging data 
reveal sub-structures in the innermost regions of 
the main disk that have never been observed be- 
fore: a 'hole' and a 'clump' are detected in both 
filters. Since these two features appear in the final 
polarized flux images (P) as well as in the frac- 
tional polarization images (pi) they are real and 
have their origin in the disk. Although these fea- 
tures are not a definite proof that planet formation 
is going on in the HD100546 disk, there is qualita- 
tive agreement between our observations and mod- 
els predicting the observational features of the in- 
teraction between low-mass planets and a disk. If 
the observed features were caused by orbiting plan- 
etary bodies in the HD 100546 disk, they should 
show a rotation of roughly 30° by 2016 (assuming 
Keplerian rotation and circular orbits). 

• Previoiis observations and models have suggested 
that a gap stretching from ~4-13 AU divides the 
circumstellar disk into an inner disk and an main 
disk. Our data support the existence of an 'inner 
rim' of the main disk around ~15 AU. At this lo- 
cation the polarized flux measured along the disk 
major axis has its maximum in both the H and Kg 
filter and it drops as we get closer to the star. 

Our data demonstrate the immense power of PDI as 
a high-contrast imaging technique at large ground-based 
telescopes to spatially resolve dusty circumstellar disks at 
small inner working angles. PDI helps to put constraints 
on the properties of the dust grains on the surface of cir- 
cumstellar disks and reveals structures within the disk 
that could arise from the interaction between the disk 
and larger bodies, such as planets, forming within the 
disk. When combined with scattered light images from 
HST (giving access to the absolute intensity of the scat- 
tered light and probing disk regions farther out due to its 
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superior surface brightness sensitivity) a comprehensive 
picture of the surface properties of circumstellar disk can 
be derived. 

In the near future, two new high-contrast instruments 
will become available both equipped with a dual-beam 
polarimeter: SPHERE at the VLT (Dohlen et al. 2006; 
Beuzit et al. 2006) and GPI at GEMINI south (Macin- 
tosh et al. 2006). Thanks to their extreme AO systems, 
these instruments will provide a much more stable PSF 
and will be able to probe polarized fiux at smaller inner 
working angles with higher sensitivity than current in- 
struments. The hope is that these instruments will mark 
a new area in resolving and analyzing the surfaces of cir- 
cumstellar disks from the ground. 

Finally, ALMA will become available as of this sum- 
mer. Once it has reached its final configuration it will 
provide vahiablc complementary information about the 
dust in the mid-plane of circumstellar disks and also 
about the gaseous component of disks with somewhat 
comparable spatial resolution. New insights - but also 
new questions - about the formation processes of planets 
lie in the not so distant future and PDI has the potential 
to play an important role. 
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Region probed by NACO/PDI data: 



HD1 00546 Main disl< surface -15-140 AU 




Inner disic Gap IVIain disic 

Ri„:-0.3AU Ri„:~13-15AU 



Fig. 1. — Sketch of the HD100546 system. Our data probe the surface layer of the main disk. 
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Fig. 2. — Final Stokes Q and U images (count rates in arbitrary units, left column) and fractional polarization pQ and pu (right column) 
of HD100546 in the H filter. The core of the PSF, where the count rate was no longer in the linear detector regime in at least one of 
the individual raw images before they were combined, has been masked out. All images have been multiplied by -1 so that ordinary and 
extraordinary beam arc the same as for the other filters to compensate for a 90° difference in the position angle of the camera (see. Table 2). 
North is up and east to the left in all images (since we aligned the disk major axis with the detector's x-axis the positive Q component is 
not aligned in North-South direction in our final images as it is the usual convention). 
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Fig. 3. — Same as Figure 2 but for the Kg filter 




Fig. 4.— Final Stokes Q images for HD100546 in the ArBl.64 filter (left) and 7B2.16 filter (right). The scaling and orientation is the 
same as in Figures 2 and 3. 
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Fig. 5. — Same as Figure 3 but for the reference star HD106797. Note that the images have the same stretch. 
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Fig. 6. — Azimuthal profiles of fractional polarization pg {H filter) and pu {Ks filter) in the top row and Stokes parameters Q {H filter) 
and U {Ks filter) in the bottom row based on Figures 2 and 3. The unit of the y-axis is '%' in the top row and 'count rate' in the bottom 
row. Curves for different annuli centered on the star are shown: diamonds (black) = 0.1-0.25", triangles (red) = 0.25—0.5", x (blue) = 
0.5—0.75", squares (green) = 0.75-1.0". Each data point shows the mean value in a 10° wedge for a given annulus. In the bottom row 
the count rates measured for the two outermost annuli have been multiplied by 3 to emphasize their sinusoidal variation. The error bars 
are the standard deviation of the mean values computed in the individual images at each dither position divided by the square root of the 
number of images that were combined (see. Table 2). 




Fig. 7. — Top: Surface brightness of polarized fiux P of the HD100546 disk in the H filter (left) and Ks filter (right). Bottom: Fractional 
polarization pj in the H filter (left) and Ks filter (right). North is up and East is left in all images. 
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Fig. 8. — Radial surface brightness profiles for the semi-major and semi-minor axes of the HD100546 disk derived from the images shown 
in Figure 7. We put a 3-pixel wide slit along the axes and computed the mean and the standard deviation across the slit for each pixel 
position. The black lines show the results for the semi-major axes, the red line those for the semi-minor axes. The blue line in the top row 
is a power-law fit to the profile of the semi-major axes as discussed in the text. The negative part of the x-axis refers to the SE semi-major 
axis and SW semi-minor axis, while the positive part shows the NW semi-major axis and NE semi-minor axis. 
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Fig. 9. — [H — ii's]-color of the polarized flux computed from the images shown in the upper row in Figure 7. North is up, east to the 
left. 
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Fig. 10. — High-pass filtered versions of the images shown in the top row of Figure 7. The images have been rotated clockwise so that 
the disk major axis runs horizontally. The white arrows indicate the positions of the two sub-structures: (1) a 'hole' and (2) a 'clump'. See 
text for more details. 
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HD1 00546 - Polarized Flux - H 
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Fig. 11. — Zoom in the inner regions of the images shown in Figure 7. The images have been rotated clockwise so that the disk major 
axis runs horizontally. The white arrow indicates the location of the "hole" (see also Figure 10). 



NACO PDI observations of the disk around HD100546 25 



100 




Fig. 12. — Resulting images from subtracting the mirrored left-hand side of the disk from the right hand side of the disk (left panel: H 
filter; right panel; Ks filter). Before subtra<;tion the images were rotated clockwise so that the disk major axis runs horizontally (same 
orientation as in Figure 11). 
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Fig. 13. — Cut along the disk major axis in the polarized flux images shown in Figure 11. Black 'x' show the results for the H images 
while red diamonds show the results for the Ks images. Pixels in the image center that were overexposed and no longer in the linear 
detector regime have been masked out. 



NACO PDI observations of the disk around HD100546 



27 



Scattering function - HDl 00546 H filter 



Scattering function - HD100546 Ks filter 
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Fig. 14. — Scattering function based on the images shown in the upper panels of Figure 7 (left: H filter; right: Ks filter). The mean 
polarized flux was computed in two annuli between 40—50 AU (black curves) and 100-110 AU (red curves) in wedges of 10° as a function of 
the scattering angle (see text). The error bars denote the standard deviation within each wedge. The flux was normalized to the maximum 
flux measured in the direction of the north-east semi-minor axis. 
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TABLE 1 
Basic parameters of objects. 



Parameter 


HD100546 


HD106797 


Reference^ 


RA (J2000) 


ll''33'"25^.44 


12''17'"06''.31 


(1) 


DEC (J2000) 


-70°11'41".24 


-65°41'34".65 


(1) 


J 


6.43±0.02 mag 


6.00±0.02 mag 


(2) 


H 


5.96±0.03 mag 


6.03±0.04 mag 


(2) 


Ks 


5.42±0.02 mag 


6.00±0.02 mag 


(2) 


Mass 


2.4±0.1 Mq 




(3) 


Age 


5...>10 Myr 


10-20 Myr 


(3),(4),(5),(6) 


Distance 


97lt PC 


96+1 


(7) 


Sp. Type 


B9Vne 


AOV 


(8) 



''References — (1) Ferryman et al. (1997), (2) 2MASS point 
source catalog (Cutri et al. 2003), (3) van den Ancker et al. (1997), 
(4) Guimaraes et al. (2006), (5) de Zeeuw et al. (1999), (6) Fuji- 
wara et al. (2009), (7) van Leeuwen (2007), (8) Houk & Cowley 
(1975). 
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TABLE 2 
Summary of observations. 



Object 


Filter 


DIT X NDIT'^ 


Dither 


Position 


Airmass 


Obs. date 


<EC)<i [%] 


(to)° [ms] 








positions'" 


angle'^ 






mean/ min /max 


mean /min/ max 


HD 100546 


H 


0.3454 s X 85 


11(15) 


155° (141.8°) 


~1.67-1.51 


2006 April 7 


36/13/47 


5.4/3.7/7.4 




Ks 


0.3454 s X 85 


12(13) 


65° (51.8°) 


~1.49-1.43 


55 


41/17/60 


4.5/3.1/7.3 




NB1.64 


3 s X 20 


4(5) 


65° (51.8°) 


~1.43 


55 


49/44/52 


4.5/4.0/5.0 




7S2.16 


0.6 s X 50 


8(9) 


65° (51.8°) 


~1.43-1.45 


55 


44/37/49 


3.4/2.6/4.3 


HD 106797 


H 


0.3454 s X 85 


1(2) 


155° (141.8°) 


~1.34-1.35 


2006 April 7 


30/23/37 


2.9/2.7/3.1 






0.3454 s X 85 


2(2) 


65° (51.8°) 


~1.33-1.34 


55 


40/40/41 


2.9/2.8/3.1 



Detector integration time (DIT) X number of integrations (NDIT), i.e., total integration time per dither position and per retarder plate position. 

^Number of dither positions used in final analysis and in parenthesis total number of observed dither position. The difference was disregarded due to poor AO 
correction. At each dither position NDIT exposures -were taken at each of the 4 different retarder plate positions (0.0°, -22.5°, -45.0°, -67.5°). 

^'Position angle of the camera on sky (north over east) as taken from the image header and after correction for the offset in the encoder value leading to a false zero 
point of the retarder plate (in parenthesis). For details see section 3.4. 

*^Average, minimum and maximum value of the coherent energy of the PSF. Calculated by the Real Time Computer of the AO system. 
^Average, minimum and maximum value of the coherence time of the atmosphere. Calculated by the Real Time Computer of the AO system. 
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TABLE 3 

HD100546 DISK INCLINATION AND POSITION ANGLE FROM FROM 

SCATTERED LKUIT IlIACER'i-. 



Reference Inclination'' PA*" Radius'' 



This work H band 


45.9° 


±1.7° 


136.0° ± 2.9° 


0.2- 


-0.5" 


This work Kg band 


48.0° 


±2.0° 


140.0° ± 2.4° 


0.2- 


-0.5" 


This work mean 


47.0° 


±2.7° 


138.0° ± 3.9° 


0.2- 


-0.5" 


Ardila et al. (2007) 


42.0° 


±5.0° 


145.0° ± 5° 


1.6- 


-2.8" 


Pantin et al. (2000) 


50.0° 


±5.0° 


127.0° ± 5° 






Grady et al. (2001) 


49.0° 


±4.0° 


127.0° ± 5° 






Augereau et al. (2001) 


51.0° 


±3.0° 


161.0° ± 5° 


0.7- 


-3.0" 



^Face-on corresponds to 0°. 

'^Position angle of disk major axis measured east of north. 
*^Used for isophot fitting. 
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TABLE Al 

Offsets used to correct for instrumental polarization in pq / pu images in the 

DIFFER EXT FILTERS. 



Object 


H 


Ks 


NB1.64 


/B2.18 


HD100546 
HD106797 


-0.01275 / -0.00049 
-0.00774 / -0.00806 


0.00463 / 0.00450 
-0.00134 / 0.00430 


-0.00768 / 0.00779 


0.01033 / 0.00995 



APPENDIX 

CORRECTION FOR INSTRUMENTAL POLARIZATION 

In the following we describe in detail how the effects of instrumental polarization were quantified and what cross- 
checks we applied to validate our approach. 

Quantification of instrumental polarization 

We use the central region around the star as zero polarization calibration source assuming that its intrinsic polar- 
ization is close to zero. An interstellar polarization component behaves like a static instrumental polarization offset 
and this would be corrected with our procedure. Assuming p « 0% for the interstellar polarization of HD 100546 
is a reasonable assumption since published polarimctry for the V-band indicate a polarization of only w 0.25 % for 
9 « 55° (Yudin & Evans 1998; Clarke ct al. 1999; Rodrigucs ct al. 2009). This value is composed of some intrinsic 
polarization component and some interstellar polarization pis{Y). Since IID100546 shows only very little extinction 
effects (Av=0.28 mag, van den Ancker et al. (1997)) we can assume that the intrinsic polarization is also low in the 
near-IR for HD100546. Of course, a real measurement would be useful but unfortunately we are not aware of any 
calibrated ncar-IR polarimctry for HD100546. 

We assumed that in the final pg and pu images the mean value of the polarization fraction averaged over an annulus 
between 4 6 pixels centered on the star is zero. In case there are no scattering dust grains around the star, and hence 
no polarization signal, the azimuthal average of the fractional polarization should be zero in any case regardless of the 
radial extend of an annulus. In case there is a disk that is not seen edge-on but with a sufficiently small inclination 
angle (which for HD100546 we know it is the case, e.g., Augereau et al. 2001) we assumed that the substructure in 
the disk is small enough so that the positive and negative contributions to the polarization fraction in the pq and 
pif images in a small centro-symmetric annulus around the star would cancel out. Ideally, if the PSF cores were not 
saturated and the central source would cancel out completely one could use the innermost pixels instead of an annulus. 
However, since our H and Ks images were saturated in the PSF core (see, section 2) we had to rely on those regions 
that were not saturated and we used the same annulus for all data sets to be consistent. In Table Al we summarize 
the mean values we found for the polarization fractions in the aforementioned annulus in the different data sets. It 
turned out that in most cases these values were small, but not zero, and they varied between all filters and both 
Stokes parameters. We assumed that these offsets were due to instrumental polarization effects and we corrected each 
individual pq and pu image for these offsets before the final images were median combined. In the final Q and U 
images these offsets arc, however, not additive offsets but factors. Hence, before we did the median combination of 
the final Q and U images we scaled the individual images with 1 + x, where x is the offset derived in the respective 
Pq and pu images. 

Validation of approach 

Once we had applied the correction factors, we checked the final results by comparing the observed structures in the 
different filters for both objects (see section 4.1). Finding strong signals for HD100546 that are similar in all filters 
and not seeing any consistent and significant structure around the reference source HD106797 gives us confidence that 
our data reduction technique is robust. 

An additional check was accomplished by analyzing the azimuthal profiles of the "butterfly" patterns (upper panels 
of Figure 6; discussed in detail in section 4.1). Comparing these figures for all flltcrs with similar figures derived 
from simple generic disk models shows that one would expect that the plotted lines computed for different annuli will 
always cross each other a,t pq = pu = irrespective of disk inclination. While this is the case in our final plots, within 
the error bars, it was not the case before the correction. 

Before- and- After comparison 

To demonstrate that instrumental polarization offset has a significant impact and needs to be corrected for we show 
in Figure Bl an uncorrected and a corrected version of the Q image of HD100546 in the H filter. As can be seen from 

Table Al the correction factor here was the largest among all images, thus this image suffered the most contamination 
by instrumental polarization. To compare these images directly to the final image shown in Figure 2 one has to 
multiply by — 1 in order to compensate for the fact that the camera has been rotated by 90° compared to the images 
in the Ks filter and hence the ordinary and extraordinary beam were switched. 
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SIGNAL-TO-NOISE ESTIMATES 

The main source of uncertainty in our images results from variations in flux in the ordinary and extraordinary beam 

between the different dither positions. This means that the Q, U, pq, and pjj images show shghtly different patterns 
in terms of signal strengths at each dither position. Compared to these variations detector read-out noise, sky noise 
and photon noise are negligible in the inner regions. Thus, it is valid to base an estimate! for the signal-to- noise (S/N) 
in our images on these variations and compare them to the final images. In addition to the median-combined final 
images, we also computed the mean images for each Stokes parameter and each filter and the corresponding sigma 
images which show the standard deviation of each pixel. The mean images and median-combined images were almost 
identical similar in terms of overall flux levels and surface brightness distribution. Since we are interested in the S/N 
of the observed extended emission and not in the S/N of individual pixels, we convolved the median images and the 
sigma images with a Gaussian kernel where the FWHM corresponded to the FWHM of the PSFs in our images. A 
S/N map can then be computed via 

S/N = VA^Dith • (Bl) 

No- 
where iVDithj I) and lo- denote the number of dither position used for the final combination, the median-combined 
image, and the sigma image, respectively. 

In Figure B2 we provide an example of these S/N maps and show those for the Q and U image in the Kg filter for 
HD100546. The orientation of the images is the same as in Figure 11 with the disk major axis running horizontally. 
To obtain the same orientation as in Figure 3 (north up and cast to the left) the images have to be rotated by 51.8° 
counterclockwise. The S/N map clearly resembles the surface brightness distribution of the disk. However, it shows 
that wherever we detect polarized scattered flux within 1.4" it is detected with a S/N>3, in particular for the total 
polarized flux P given by P = y/Q-^ + C/2. 

For comparison, we show the S/N map of the U image in the Kg filter for HD106797 in Figure B3. Note the different 
image stretch compared to the S/N maps of HD100546. Here, no feature is detected with a S/N>1.5. 
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Fig. B2. — S/N maps of the Q (left) and U (right) images for HD100546 obtained in the Kg filter. The disk major axis of the disk runs 
horizontally. The color scale runs from to 12. See text for more details. 
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Fig. B3. — S/N map of the U image for HD106797 in the Ks filter. Note that the color scale runs from to roughly 1.5. 



